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that  is  equivalent  to  the  second  order  stochastic  dominance  relation 
when  certain  independence  conditions  hold.  An  example  illustrates  these 
results. 
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ABSTRACT:  In  this  paper,  we  examine  a  single-stage,  multiobjec¬ 

tive  decisionmaking  problem  under  uncertainty.  The  decisionmaker 
can  select  any  one  of  a  finite  number  of  alternatives.  After  any 
alternative  is  chosen,  one  of  a  finite  number  of  outcomes  will 
result.  The  probabilistic  relationship  between  each  alternative 
and  each  outcome  is  presumed  to  be  known.  We  assume  that  all 
that  is  known  about  the  decisionmaker  is  that  he  or  she  is  risk 
averse.  Our  objective  is  to  determine  the  smallest  subset  of 
alternatives  that  is  guaranteed  to  contain  the  most  preferred 
alternative  on  the  basis  of  this  assumption.  The  achievement 
of  this  objective  presumably  enhances  decisionmaking  since  alter¬ 
native  selection  is  generally  easier  if  made  from  a  subset  of 
the  alternative  set  rather  than  from  the  entire  alternative  set. 

The  intent  of  this  paper  is  to  present  an  approach  which 
achieves  this  objective  and  which  has  computational  times 
amenable  to  interactive  decision  aiding.  We  make  use  of  a  fact, 
due  to  Fishburn  and  Vickson, which  states  that  the  feasibility 
of  a  certain  collection  of  linear  equalities  and  inequalities 
represents  a  necessary  and  sufficient  condition  for  one  alterna¬ 
tive  to  be  weakly  preferred  to  another  with  respect  to  the  second 
order  stochastic  dominance  (SSD)  relation.  The  approach  presented 
here  uses  transitivity  and  upper  and  lower  bounds  on  this  rela¬ 
tion  in  order  to  reduce  the  number  of  concomitant  linear  programs 
necessary  for  solution.  The  lower  bound  is  provided  by  the  first 
order  stochastic  dominance  relation;  the  upper  bound  is  given  by 
a  relation  that  is  equivalent  to  the  second  order  stochastic 
dominance  relation  when  certain  independence  conditions  hold.  An 
example  illustrates  these  results. 


♦This  research  has  been  supported  by  ONR  Contract  N0014-80-C-0542 . 


INTRODUCTION 


A  well-known  approach  for  multiobjective,  single-stage 
decision  aiding  under  uncertainty,  c.f.  (Keeney  and  Raiffa, 

1976),  has  the  following  tasks  associated  with  it: 

1)  Determine  the  utility  function  of  the  decisionmaker  (DM1 

2)  Calculate  the  expected  utility  for  each  alternative. 

(We  assume  throughout  that  the  outcome  probabilities 
as  a  function  of  alternative  are  given.) 

3)  Select  the  alternative  having  the  largest  expected 
utility. 

A  major  difficulty  with  implementing  this  approach  in  practice 
is  that  utility  function  assessment  is  often  a  stressful 
task  requiring  a  substantial  amount  of  time  and  effort.  Addi¬ 
tionally,  utility  assessment  can  require  cognitive  perspectives 
not  within  the  previous  experience  of  the  DM,  which  can  pro¬ 
duce  results  of  lower  quality  than  potentially  achievable.  A 
potentially  useful  tactic  for  reducing  these  difficulties  is 
to  investigate  the  implications  of  a  less  than  complete  descrip¬ 
tion  of  the  utility  function,  e.g.  the  various  stochastic  domi¬ 
nance  procedures  (Fishburn  and  Vickson,  1978) . 

Less  than  a  complete  description  of  the  utility  function, 
however,  almost  invariably  produces  less  than  a  total  ordering 
on  the  alternative  set,  and  a  weaker,  partial  ordering  on  the 
alternative  set  usually  cannot  identify  the  most  preferred  alter- 
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native.  A  partial  order  can,  however,  identify  the  nondom- 
inated  set,  a  set  of  alternatives  that  is  guaranteed  under 
mild  conditions  to  contain  the  most  preferred  alternative.  The 
identification  of  the  nondominated  set  is  often  quite 
adequate  for  decision  aiding;  e.g.  see  (White  et.  al,  1980) 
for  a  medical  decisionmaking  situation.  We  remark  that  the 
number  of  alternatives  in  the  nondominated  set  is  dependent  on 
the  amount  of  partial  preference  information  known;  see  (White 
and  Sage,  1980a,  1980b)  for  further  discussion. 

The  apparent  fact  that  the  nondominated  set  is  often  a 
sufficiently  informative  aid  for  decisionmaking  has  motivated 
the  development  of  decision  aiding  procedures  that  allow  the 
mix  of  alternative  order  specificity  and  utility  (or  value) 
function  identification  time,  effort  and  stress  to  be 
adaptively  determined  by  the  DM  (White  and  Sage,  1980a,  1980b). 
Adaptive  determination  of  this  mix  requires  that  the  time 
necessary  to  determine  the  impact  of  additional  (or  initial) 
preference  information  on  alternative  order  specificity  be  small 
due  to  the  substantial  constraints  that  are  often  placed  on  the 
available  time  of  most  DM's.  One  of  the  decision  aiding  pro¬ 
cedures  presented  in  (White  and  Sage,  1980b)  is  based  on  second 
order  stochastic  dominance  (SSD)  and  hence  on  only  the  often 
quite  behaviorally  relevant  assumption  that  the  DM  is  risk 
averse.  A  straightforward  implementation  of  this  SSD-based 
decision  aiding  procedure  requires  the  formulation  and  solution 
of  P(P-l)  linear  programs  for  a  problem  having  P  alternatives. 
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Our  current  computational  experience  indicates  that  such  an 
implementation  almost  invariably  requires  computation  times 
that  are  unacceptably  large,  even  for  small  problems,  for  inter¬ 
active  decision  aiding.  The  intent  of  this  paper  is  to 
propose  a  set  of  procedures  that  may  significantly  reduce 
the  computational  time  and  effort  necessary  to  order  the 
alternative  set  using  SSD,  thus  enhancing  the  potential  of 
the  SSD-based  approach  as  an  interactive  decision  aiding  pro¬ 
cedure. 

This  paper  is  organized  as  follows.  The  problem  is 
formulated  and  preliminary  results  and  definitions  are  given 
in  Section  2.  In  Section  3,  we  state  and  investigate  condi¬ 
tions  that  may  often  reduce  the  time  required  by  the  approach 
taken  in  (White  and  Sage,  1980b)  to  order  the  alternative 
set  using  SSD.  Conclusions  are  presented  in  the  final  section. 


PROBLEM  FORMULATION  AND  PRELIMINARIES 


We  assume  that  the  DM  can  select  for  implementation  any 
one  of  P  predetermined  alternatives  from  the  alternative  set 
n=  {it  ,..,it  }.  After  an  alternative  is  implemented,  any  one 
of  M  possible  outcomes  will  occur.  There  are  N  objectives 
under  consideration.  Let  v™  be  the  predetermined  value  score 
of  the  in  outcome  with  respect  to  the  n1"  objective.  The 
real  number  v™  is  isotone  (monotonically  nondecreasing)  in 
preference  with  respect  to  the  n1"  objective;  that  is,  out¬ 
come  m'  is  (weakly)  preferred  to  outcome  m  with  respect  to 
objective  n  if  and  only  if  vm  >  vm.  Let  V  =  {vm,  m  =  1,...,  M), 
the  set  of  all  value  score  vectors,  where  vm  =  {v™,...,  v™}. 

The  probability  that  outcome  m  will  result  if  alternative  irP 
is  selected  is  irp(vm)  . 

We  assume  that  there  exists  a  (presumably  unassessed) 
utility  function  u:  V  R  which  reflects  the  DM's  preferences 
in  that  outcome  m'  is  (weakly)  preferred  to  outcome  m  with  all 
objectives  under  consideration  if  and  only  if  u(v  )  >  u(vm). 
Alternatives  are  compared  on  the  basis  of  expected  utility: 
alternative  ir '  is  (weakly)  preferred  to  alternative  it  if  and 
only  if  E(u,  tt  ' )  >_  E(u,  it)  ,  where 

E  (u,  it)  =  2  u  (v)  tt  (v)  . 

veV 

We  make  the  following  two  assumptions  about  the  DM: 


4 


1)  Assume  that  the  DM  prefers  outcome  m'  to  outcome  m 
when  only  objective  n  is  considered.  Let  this 
statement  be  true  for  all  n  =  1,..,  N.  Then,  the 
DM  prefers  outcome  m'  to  outcome  m  when  considering 
all  objectives  simultaneously  (consistency) . 

2)  The  DM  prefers  the  expected  consequence  of  a  lottery 
(the  concept  of  a  lottery  is  discussed  at  length  in 
(Keeney  and  Raiffa,  1976))  to  that  lottery  (risk 
aversion) . 

These  two  assumptions  imply  that  the  DM's  utility  function  is 
isotone  and  concave,  respectively,  and  thus  is  a  member  of  the 
set  U2  =  {u:  u  is  isotone  and  concave).  The  objective  is  to 
provide  the  DM  with  a  set  of  alternatives  which  is  guaranteed 
to  contain  the  most  preferred  alternative  under  the  assumption 
that  all  that  is  known  about  the  DM's  utility  function  is  that 
it  is  a  member  of  l^. 

We  say  that  alternative  ir '  is  (weakly)  preferred  to  alter¬ 
native  it  with  respect  to  SSD,  i.e.  if  and  only  if  E(u,n')>^ 

E(u,  it)  for  all  u  e  .  Thus,  if  the  DM  is  consistent  and  risk 
averse,  tt  '  R2  w  implies  that  w '  can  be  expected  to  be  at  least 
as  good  an  alternative  choice  as  it.  It  is  well-known,  c.f. 

(White  and  Sage,  1980a) ,  that  the  most  preferred  alternative  is 
a  member  of  the  nondominated  set  (it  e  II:  there  does  not  exist 
a  tt  '  e  n  such  that  tt'K2  it  and  not  n  R 2 tt  '  > .  Determinincr  the 
nondominated  set  of  n  with  respect  to  R2,  or  more  generally 
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determining  the  set  P2Q  n  *  where  (it',  tt )  e  ?2  if  and  only  if 
iv'  l?2  w,  helps  to  restrict  the  search  for  the  most  preferred 
alternative,  thus  presumably  aiding  decisionmaking.  (Note  that 
ir  is  nondominated  if  there  is  no  i1  e  n  such  that  (ir',ir)e  P2 
and  (tt,  tt ')/  P2;  thus,  knowledge  of  P2  can  used  to  determine 
the  nondominated  set  of  IT  with  respect  to  P2>) 

The  following  necessary  and  sufficient  conditions,  a 
slightly  generalized  version  of  a  result  presented  in  (Fishburn 
and  Vickson,  1978) ,  suggest  an  approach  for  determining  P2« 

THEOREM  1:  it*  R2  it  if  and  only  if  there  exists  a  feasible 
solution  to  the  set  of  linear  equalities  and  inequalities: 

(i)  d„  >_  0  for  all  i,  j  =  1,...,  M 

M 

(ii)  Z  d..  =  1  for  all  i  =  1,...,  M  such  that  ir’fv1)/  0 

j  =  l  13 

M 

(iii)  ^(v-5)  =  Z  ir'iv1)  d..  for  all  j  =  1,...,  M 
i=l  1J 


(iv) 


M 

Z 

i=l 


d. ,  v^  <  v1  for  all  i 
lj  n  —  n 


=  1,. 


M  such  that 


it  '  (v1)  /  0  and  all  n  =  1,..., 


N. 


A  probabilistic  interpretation  of  {d^ }  is  given  in 
(Fishburn  and  Vickson,  1978)  . 

A  brute-force  application  of  Theorem  1  for  determining  P2 
requires  the  formulation  and  solution  of  P(P-l)  linear  programs. 


-  9 
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each  having  up  to  M (M  +  N  +  2)  decision  variables,  2M  of  which 
are  artificial  variables,  and  up  to  M (N  +  2)  side  constraints. 
Our  present  level  of  experience  with  the  process  of  formu¬ 
lating  and  solving  these  linear  programs  indicates  that  they 
may  often  require  computer  time  that  is  unacceptably  large 
for  interactive  decision  aiding,  even  with  relatively  small 
problems  and  a  relatively  fast  computer.  For  example,  the 
P=6,  M=5,  and  N=4  problem  considered  in  Example  3 (White  and 
Sage,  1980b)  required  roughly  1.5  minutes  of  CPU  time  on  the 
CDC-6400  at  the  Computing  Center  of  the  University  of  Virginia. 
Since  our  computations  were  done  interactively  in  a  time 
sharing  mode;  turn  around  time  can  be  expected,  and  was 
experienced,  to  be  between  5  and  15  minutes,  depending  on  the 
system  load.  This  length  of  time  seems  excessive  for  inter¬ 
active  decision  aiding.  In  the  next  section,  we  present  tech¬ 
niques  that  may  often  reduce,  sometimes  substantially,  this 
computational  burden. 

MAIN  RESULTS 

In  the  previous  section,  a  procedure  for  determining  P2 
was  suggested  that  involved  a  straightforward  application  of 
Theorem  1.  In  this  section,  we  present  several  results  that 
may  often  reduce,  sometimes  substantially,  the  computational 
times  associated  with  this  procedure.  These  results  are 
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presented  and  proved  following  three  preliminary  defini¬ 
tions. 

1.  The  alternative  it'  is  said  to  be  (weakly)  preferred 
to  ir  with  respect  to  first  order  stochastic  dominance  (FSD) , 
i.e.  ir'R^  tt,  if  and  only  if  E(u,  it')  >_  E(u,  it)  for  all  ueU^= 

{u:  u  is  isotone}. 

2.  The  alternative  tt  *  is  said  to  be  (weakly)  preferred 
to  tt  with  respect  to  strong-SSD  (SSD)  ,  i.e.  tt'J^  11  >  if  and 
only  if,  for  each  n  =  1,...,  N,  there  exists  a  feasible  solution 
to  the  set  of  linear  equalities  and  inequalities  (i) ,  (ii) , 

(iii) ,  and 

M 

(iv)  '  £  d .  ,  v-1  <  v  for  all  i=l,...,  M  such 

i=l  ij  n  -  n 

that  -rr '  (v1)^  0 . 

3.  An  arbitrary  relation  R^  on  IT  is  said  to  be  stronger 
than  (more  precisely,  at  least  as  strong  as)  an  arbitrary 
relation  R  on  IT,  i.e.  R  C  R,  ,  if  and  only  if  it  '  R  tt  implies 
tt  '  Rfo  tt,  for  any  pair  tt  '  ,  tt  e  n. 

We  remark  that  tt'R2  tt  is  equivalent  to  N  separate  checks 
for  univariate  SSD,  for  which  there  exists  a  computationally 
simple  procedure,  c.f.  Section  2.14  in  (Fishburn  and  Vickson, 
1978).  We  also  note  that  under  certain  independence  conditions 
(presented,  for  example,  in  Theorem  2 . 11  of  (Fishburn  and  Vickson, 
1978)),  R^  =  R" 2 ,  i.e.  R2  C  R2  and  1<2  Q.  f^*  We  noW  Present  our 
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main  result,  the  proof  of  which  is  in  the  Appendix. 


THEOREM  2:  R^,  R2,and  are  transitive,  and  R^C  R2  Q.^2' 

The  impact  of  this  result  is  due  to  the  fact  that  it 
can  be  used  to  reduce,  often  drastically,  the  number  of 
linear  programs  that  require  solution  in  order  to  construct 
the  set  ?2 •  Transitivity  implies  that  it  is  not  necessary 
to  check  whether  or  not  it"  R  it  ,  if  n"  R  tt  *  and  -n '  R  it,  for 
arbitrary  transitive  relation  R.  R^  C  R^  C  R^  implies  that 

once  P^  and  P2  are  known  (where  we  define  P^  and  P2  similarly 
to  P2) f  the  only  pairs  that  require  the  linear  program  check 
in  order  to  complete  determination  of  P2  are  those  pairs  in 
?2  which  are  not  in  P^  and  which  have  not  already  been  added 
to  ?2  by  the  above  transitivity  argument.  Theorem  2,  there¬ 
fore,  suggests  the  following  four  step  procedure  for  determining 


(1)  Determine  P^.  (Relatively  computationally  simple 
procedures  for  determining  P-^  are  presented  in  (White 
and  Sage,  1980b).) 

(2)  Determine  P.,.  (Relatively  computationally  simple  pro¬ 
cedures  for  determining  P2  are  suggested  in  (Fishburn 
and  Vickson,  1978).) 

(3)  Evaluate  all  pairs  in  P  which  are  not  in  using 

Theorem  1. 


(4)  Construct  P2  by  adding  the  appropriate  pairs  from  Step  3 
to 


The  above  procedure  prompts  three  comments.  First, 
the  transitivity  of  all  relations  can  be  useful  in  reducing 
the  number  of  linear  program  formulations  and  solutions, 
a  fact  not  explicitly  mentioned  above.  Second,  it  is  incon¬ 
sequential  that  Step  1  is  performed  prior  to  Step  2.  Which¬ 
ever  of  the  first  two  steps  is  performed  first  can,  however, 
impact  on  the  time  required  to  perform  the  second  step.  This 
impact  is  due  to  the  facts  that  if  (it*,  tt)  e  P^,  then  ( -it •  ,  n)  e 
?2  and  if  (tt',  n)  /  P.^ ,  then  (it1,  it)  /  P^,  since  P^Q^P^. 

Third,  in  checking  for  strong-SSD,  it  is  necessary  to  examine 
only  down  to  the  first  objective  that  fails  to  satisfy  the 
univariate  SSD  criterion  (if  one  exists) . 

We  now  present  an  example  illustrating  the  above  proce¬ 
dure. 

EXAMPLE:  Consider  Example  3  in  (White  and  Sage,  1980b).  In 

that  example  problem,  there  were  six  available  alternatives, 
five  possible  outcomes,  and  four  objectives  under  consideration 
prior  to  the  objective  aggregation  procedure,  i.e.  P=6,  M=5,  and 
N=4 .  The  Table  presents  the  assumed  data.  Results  in  (White 
and  Sage,  1980b)  indicate  that  P^  =  {(4,3)}.  Calculations  based 
on  procedures  suggested  in  Section  2. 14  of(Fishburn  and  Vickson, 
1978)  show  that  T ^  =  ((1,3),  (2,3),  (4,3)}.  Thus,  it  is  only 

necessary  to  check  the  pairs  (1,  3)  and  (2,  3)  in  order  to 
determine  P2-  Solution  of  the  two  associated  linear  programs 


Outcome  Number 


Objective 

Number 


(a)  Value  Scores  for  Each  Outcome  and  Objective 


Outcome  Number 


3 

II 

<N 

II 

e 

m=3 

m=4 

m=5 

P=1 

0-6 

0.1 

0.2 

0.1 

0.0 

p=2 

0.7 

0.0 

0.1 

0.2 

0.0 

p=3 

0.3 

0.1 

0.0 

0.4 

0.2 

p=4 

0.3 

0.0 

0.1 

0.1 

0.5 

p=5 

0.1 

0.1 

0.0 

0.1 

0.7 

p=6 

0.0 

0.1 

0.1 

0.0 

0.8 

(b)  Outcome  Probabilities  for  Each  Alternative 


Table: 


Data  for  the  Example 


indicates  that  (1,  3)  e  ?2  and  (2,  3)  i  ? 2>  and  hence  ?2= 

P1t/  (1,  3)  =  {(1,  3),  (4,  3)}.  Note  that  P2*V2m  This 
result  is  in  agreement  with  a  result  found  in  (White  and 
Sage,  1978) ,  which  was  determined  from  the  formulation  and 
solution  of  P(P-l)  =>  6(6-1)  =  30  linear  programs. 

The  first  three  objectives  were  linearly  aggregated  in 
Example  3 (White  and  Sage,  1980b)  with  weights  0.1,  0.1, 
and  0.8,  respectively.  As  a  result,  P^  =  {(1,2),  (4,3), 
(6,5)}.  Calculations  show  that  ~P2=  {(1,2),  (1,3),  (2,3), 

(4.3) ,  (6,5),  (6,3),  (5,3)}  (note  that  (1,3)  and  (6,3)  are 

members  of  F2  by  transitivity) .  Thus,  the  only  pairs  that 
require  examination  by  the  procedure  suggested  in  Theorem  1 
are:  (1,3),  (2,3),  (6,3),  (5,3).  We  observe  that  if  (2,3)e 

?2  and  (5,3)  e  P2 ,  then  it  is  not  necessary  to  check  if 

(1.3)  e  ?2  and  (6,3)  e  P2,  respectively,  because  of  the  trans 
itivity  of  P2.  Solution  of  the  associated  linear  programs 
show  that  (1,3)  and  (6,3)  are  members  of  P2.  We  have  there¬ 
fore  determined  that  P 2 =  V 2  by  formulating  and  solving  only 
two  linear  programs.  This  result  is  in  agreement  with  a 
result  found  in  (White  and  Sage,  1980b) ,  which  again  was  deter 
mined  from  the  formulation  and  solution  of  30  linear  programs 


CONCLUSIONS 


This  paper  has  investigated  procedures  for  making 
SSD  a  viable  concept  for  interactive  decision  aiding.  Our 
primary  contribution  toward  achieving  this  objective  has 
been  the  identification  of  a  partial  order  that  acts  as  an 
upper  bound  on  the  SSD  partial  order.  Our  present  level 
of  experience  indicates  that  this  upper  bound,  the  FSD 
partial  order  lower  bound,  the  transitivity  of  all  three 
partial  orders,  and  necessary  and  sufficient  conditions 
due  to  (Fishburn  and  Vickson,  1978)  can  often  be  used  to 
obtain  a  significant  reduction  in  the  computational  demands 
associated  with  a  straightforward  application  of  Theorem  1 
in  determining  P2- 

The  four  step  procedure  for  determining  P 2  proposed 
here,  however,  may  not  always  reduce  computational  time. 

Although  straightforward  application  of  Theorem  1  requires 
formulating  and  computing  at  least  as  many  linear  programs 
as  required  by  the  four  step  procedure;  it  does  not  require 
the  determination  of  P^,  or  transitivity  checks.  If 

P^  =  >)>  and  P2  =  n  x  n,  it  is  clear  that  an  application  of 

Theorem  1  that  allows  for  transitivity  checks  will  be  computa¬ 
tionally  quicker  and  therefore  superior  to  the  four  step  pro¬ 
cedure.  We  have  found,  however,  that  the  number  of  pairs  in  wh 
are  not  in  P^  have  usmlly  been  a  ana  11  fraction  of  P(P-l)  and  that  the  time 


necessary  to  calculate  and  T ^  has  typically  been  signi¬ 

ficantly  smaller  than  the  time  necessary  to  formulate  and 
calculate  the  additional  linear  programs.  Future  computa¬ 


tional  experience  is  expected  to  further  indicate  the  merits 
of  both  approaches  for  calculating 
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Appendix:  Proof  of  Theorem  2. 

It  is  shown  in  (Fishburn  and  Vickson,  1978)  that 
R^  and  R2  are  transitive  and  that  R^ £  K2.  A  simple  argu¬ 
ment,  based  on  the  fact  the  univariate  SSD  relation  is 
transitive,  proves  that  f?2  is  transitive.  In  order  to  prove 
R2 Q  %2>  define  VQ  as  the  set  of  all  { ^ }  such  that  (i) , 

(ii)  ,  and  (iii)  hold  and  as  the  set  of  all  {d^  }  such  that 
(iv) '  holds.  Note  that: 

(a)  it1  R2  it  if  and  only  if 


(b)  tt*  R2  tt  if  and  only  if 

VQC\  t>n  /  <f>  for  all  n  =  1,...,  N. 
Use  of  the  fact  that 

’."[A’.]- 

easily  implies  K2~  R2* 
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